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I. INTRODUCTION 


The detection and location of submarines (and other 
magnetic bodies) through the discrimination of changes. or 
anomalies in the Earth's magnetic field is called Magnetic 
Anomaly Detection or MAD. In this technique, a magnetometer 
measures the magnitude of the Earth's magnetic field and 
provides an indication of that magnitude, or, more usually, 
an indication of changes in the magnitude of the Earth's 
field. These changes, or anomalies, can indicate the 
presence of magnetized bodies which may or may not be a 
submarine, 

The magnetometer currently in use in the United States 
Navy for use in this MAD process is the AN/ASQ-81 
metastable helium vapor total field magnetometer. 

Research is currently being conducted at the Naval 
Postgraduate School in Monterey, California, in various 
aspects of the applications of magnetometers, including 
Magnetic Anomaly Detection (MAD). Within the context of 
this research, magnetic field measurements are made through 
the use of sets of wire wound coils vice any specific 
Magnetometer or magnetic detecting system. The data 


collected through the use of these coils is evaluated and 
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processed in a variety of methods for different project 
goals. 

This thesis project is designed to produce an acceptable 
alternative to the physical presence of an experimental 
AN/ASQ-81 Magnetometer at the postgraduate school by 
allowing the determination, in conjunction with ether 
research in progress, of the output of the AN/ASQ-81 
magnetometer from the data collected from the_ school's 
measurement coils. It is hoped that this will assist future 
research projects asi, for example, in allowing a 
determination of environmental noise of such characteristics 
as to affect the AN/ASQ-81 magnetometer operationally with 
the eventual goal of providing an environmental noise index 
or a system of removing such noise from the magnetometer- 


detection system. 
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II. GEOMAGNETICS REVIEW 


A. EARTH'S MAGNETIC FIELD 


1. Constituents of the Geomagnetic Field 

The most common method of specifying the consituting 
parts of the geomagnetic field is to divide the field in 
terms of distance from the center of the Earth. This 
method results abil three clasisifications: internal, 
Grustal, and external. [Ref. 1] 

The internal field originates in the core region and 
is the most stable field, containing only extremely low 
frequency temporal variations. The crustal, or anomalous, 
field arises from modifications made on the internal field 
Meematerials and structures in the Earth's crust. These 
variations are not constant with regard to spatial 
locations, and comprise part of what is known as geological 
variations. The external field is the most dynamic and 
arises from many sources, including the interaction between 
the solar wind and the Earth's magnetic field. 

In addition to this method of defining the Earth's 
Magnetic field is the method of time variations. This 
method consists of considering that part of the field which 


varies with periodicities greater than about one year as the 
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steady field and everything else as the variation field. 
(eet. 2] 

The steady field consists of the internal field, also 
refered to as the main field. Slow variations of the main 
field with periods of years or longer are referred to as 
secular variations. 

There are various elements that contribute to the 
geomagnetic field, some of which are external to the Earth's 
surface. External contributions make up only a small part 
of the steady field, but play an important role in the 
variation field. These external sources include current 
systems in the Earth's upper atmosphere affected by _ solar 
electromagnetic radiation and gravitation, solar corpuscular 
radiation and the interaction of solar plasma with the main 
Preid, and the effect of the solar inteplanetary field. 
[Ref. 3] 

The geomagnetic field changes with time. As 
previously mentioned, very slow variations in the main field 
with periods of on the order of years to thousands of years 
are referred to as secular variations. Secular variations 
are caused by a variation in the strength or orientation of 
the Earth's center dipole. 

Other time variations of the field can be categorized 


into quiet variation fields and disturbed variation fields. 


Disturbed variation fields include geomagnetic 
micropulsations, which are of particular interest to 


12 


“a 
AEN... 
Blei> pajieszey S43 980 ARES Y VOTe ae 


we 
= =" 
o 

; 


: : ir... 
beis jie)? Leneeast ods Se leeeseneS elel? yhees 


lh. avn? eivey volz a3 ofes : oo ‘ 
_ : : ° oe 
” -s 2Te cagtel 5 eveme to ebotveq. “ag 
, 7 . 


_ ; 
itseitae 


cA | 


ed 


ud.cerpe@s ted? eafteegie deedeany . a0 


aa. oe 
: ‘2 99 [etastys 478 SQPGm 2@ Sade +bleit 923 
noc e ‘len qv c@am an@gateJasneo ~Lantesed; @& 
7 eee Y 
ass qe we hie Ged ,biels cose |. 
a 
cers setuda eminot danyveéeee eesdT ihletag 


a, “fe ; 
‘ir stedqteere 6 seed © 'ASseD add 2 
*-# eos 

itt pcfeats.tve soebinthand st 


mai 
trice te onfaseyean offs bas lm 


oe 


of? ban 
——— . 


di 


LOE 3 : 72atis 


t ss 


‘> sfeangseeey 6 lett 


: BV ols . beterraen 21 
bad issy°% Teht° en? @p Fo nook 


iW3ee @8 67 teste 
stevie ed¢ ol selesiaes o ed beeen 
2 pqtbh Gatas> e'dsamq 


] 

6634 ef7 l¢ enessarra¥’ GOES Gadst. 

=o Ui -:fites 280g teeh San shies? geese eee eer 
on Shia z)— ett » (lanl way ‘ bedt 


i465. Io o78@—e5a8e 


operational forces as these can mask target signatures) and 
are therefore a source of noise to MAD sensors. 

Quiet variation fields are those which are not due to 
disturbances in the interplanetary environment and which 
vary slowly and regularly. [Ref. 3] 

Disturbed variation fields are geomagnetic field 


variations that appear to be the result of interplanetary 


environmental changes and do \not posess a simple 
periodicity. These variations include ionospheric 
disturbances, the aurora, geomagnetic storms, and 


geomagnetic micropulsations. 
2. Elements of the Magnetic Field Vector 

The geomagnetic field vector is characterized at any 
point by its direction and magnitude. This is commonly 
accomplished through a system of coordinates as shown in 
Pagure 2.1. The field is measured in terms of local 
coordinates with respect to true North. [Ref. 3] 

The various coordinates are referred to as magnetic 
elements and are defined as follows: 

B: Total field intensity (the symbol F is 


sometimes also used, as in this figure.) 


H: Horizontal component 

X: Northward, or NorthSouth component 
Y: Eastward, or EastWest component 

Z: Downward, or Vertical component 
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Figure 2.1 Magnetic Field Elements [Ref. 4]. 


D: Declination or magnetic variation 


This is the angle between X and H and 


is measured positive eastward. 


i tnelination or Dip Angle. 


This is the angle between H and B (or F) 


and is measured positive downward. 
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TIT. THE AN/ASQ-81 MAGNETOMETER 


A. DESCRIPTION OF SYSTEM OPERATION 

The Magnetic Anomaly Detecting set currently in use in 
the U S Navy is the AN/ASQ-81 magnetometer. This set is 
used to locate and classify submerged submarines by sensing 
disturbances in the Earth's magnetic field (anomalies) 
caused by the presence of the magnetic mass of the 
submarine. The disturbance of the Earth's field is detected 
by the magnetometer, processed through filtering circuits, 
and amplified. The output signal of the magnetometer is 
@asplayed on a chart recorder for interpretation by an 
operator. 

The magnetic detecting set is a metastable helium vapor 
Magnetometer. The operation of the magnetometer is based on 
the light absorbtion properties of helium gas subjected to 
certain light stimulus (optical pumping), radio frequency 
excitation, and the Earth's magnetic fied. The 
Magnetometer consists of a helium lamp, lens and polarizer 
to generate a beam of polarized light radiation. This 
focused and polarized light beam is directed through a 
Heitium absorbtion cell to an infrared (IR) detector. Some 
of the helium gas in the absorbtion cell is maintained in a 


metastable state by application of VHF excitation. 
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Figure 3.1 : Metastable Helium Magnetometer [Ref,5]. 


The Earth's magnetic field imposes a magnetic force upon 
the excited helium vapor atoms to force the atoms into one 
of three energy sublevels. This is called the Zeeman 
effect. The rate or frequency of atomic precession caused 
by this effect is called the Larmor frequency. A helium 
lamp is used to optically pump the atoms in the absorbtion 
Cell, with the result that the- polarized light energy 
passing through the absorbtion cell will polarize 
(magnetize) the helium atoms in the absorbtion cells. by 
selectively pumping the Zeeman levels of the energy of the 
helium atoms in the cell. The magnetization direction is 
determined by the polarization of the photons from the 


helium lamp. 
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RF energy is then introduced to the absorbtion cell in 
the form of an additional magnetic field imposed through the 
use of coils oriented perpendicular to the precessed 
polarized helium atoms in the absorbtion cell and energized 


by a variable frequency RF oscillator. The RF oscillator is 


tuned to the Larmor frequency, which results in 
depolarization of the atoms. The atoms attempt to equally 
repopulate the Zeeman energy levels. However, the helium 


Hamp is still beaming polarized light energy into the 
absorbtion cell, causing the atoms to absorb light energy 
and rise to an excited energy levbel. This absorbtion of 
light energy is detected through the use of an infrared 
detector. The RF oscillator frequency producing maximum 
light absorbtion is called the resonant frequency, and is 
determined through the use of a servo loop from the infrared 
detector to the RF variable frequency oscillator. 
Therefore, any change in the Earth's magnetic field 
intensity will result in a change in the Larmor frequency of 
the helium atoms in the helium absorbtion cell. This new 
Larmor frequency will be detected by the ASQ-81 
magnetometer. Since the gyromagnetic ratio of helium is 
28.024 HZ per gamma, this detection of the resonant 
frequency provides a measurement of the Earth's magnetic 


field intensity at any given time. A change in the Earth's 
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magnetic field intensity could signal the presence of a 
submerged submarine. 

The output resonant frequency developed by the 
Magnetometer is converted to a proportional output voltage 
which is filtered through the Magnetic Anomaly Detection 
(MAD) bandpass filters for environmental noise reduction 
and utilized to drive a chart recorder for observation by 


an operator. [Ref. 6] 


B. TRANSFER FUNCTIONS 

Transfer functions for the AN/ASQ-81 filters were 
obtained from the manufacturer of the AN/ASQ-81 detecting 
set, Texas Instruments of Dallas, Texas. These transfer 


functions are listed in Appendix A and are in the form of 


H(s), that is, the frequency domain, or S domain, where 
S = jw. The s-domain representation for transfer functions 
is routinely utilized to express output system 
characteristics for given system inputs. As the S domain 


representation is not utilzed further in this discussion, it 
will not be further explained. 

As the output signal of the ASQ-81 magnetometer is 
filtered through a fixed high-pass system, then through a 
selectable low pass system and a selectable high pass system 
(as shown in Figure 3.2 below), the transfer funtions are 


listed in this order. 
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Further discussion will be made of the selectable filters 


later. 
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IV. DIGITAL FILTERING MODELLING OF SYSTEMS 


A. SEQUENCE REPRESENTATION OF TIME FUNCTIONS 
1. Signal Representation 

A signal can be defined as a function that conveys 
information, generally about the state or behavior of a 
physical system. Although signals can be represented in 
many ways, the information conveyed by the signal is 
contained in a pattern of variations of some form. Signals 
are represented mathematically as functions of one or’ more 
independent variables, one of the most common of which is 
time. 

The independent variable of the mathematical 
representation of a signal may be continuous or discrete. 
Continuous time signals are signals that are defined over 
continualy values of time and are therefore represented by 
continuous-variabled functions. Discrete time signals are 
defined at discrete time intervals and are therefore 
represented by functions whose independent variable(s) take 
on discrete values only. Discrete-time signals are 


represented as sequences of numbers. [Ref. 7] 


In addition to ?the “fact\ that the independent 


variables can be either continuous or discrete, the signal 
amplitude can be either continuous or discrete. Digital 
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signals are those for which both time and amplitude are 
discrete. Analog signals are those for which both time and 
amplitude are continuous. 

Digital signal processing deals with transformations 
of signals that are discrete in both time and amplitude, 
usually represented by sequences of numbers. The nth number 
in the sequence x being processed is usually represented as 
mon), and is formally written as: 

x=[x(n)], -Oo< n < +€ 


In general, an arbitrary sequence can be expressed as 


CO 
eC) = x(k) dC(na=k) 
k= -@& 
where d(n-k) is the unit sample at time k. In other words, 


an arbitrary sequence may be expressed as a sum of scaled, 
shifted unit samples, where the scaling factor is equal to 
the amplitude of the sequence at that time. 
2. Linear Shift-Invariant Systems 

A system is defined mathematically as a unique 
transformation or operator that maps an input sequence 
[x(n)] into an output sequence [y(n)]. This is denoted 
as: 

y(a)\= T[Lx(n) ] 


aad is often depicted as in Figure 4.1. 
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Figure 4.1: Representation of Transformation of an Input 
Sequence to an Output Sequence. [Ref.7] 
Classes of discrete time systems are defined by 
placing constraints on the transformation T[]. 
thesfclassy,of,.Jinear, systems -ais. defined by the 
Brancaple of superposition. Te we( ny), and. y. (n)4.are;y the 
il 2 
responses when x (n) and x (n) are the inputs, then a system 
ia linear if 
Dhax Gu) 4. bx Cn) j= allx Cn)] + bilx @)] 
1 Zz 1 2 
= aye (a) by. (a> 
if 2 
for any arbitrary constants a and b. This,. together with 
the concept of representing a sequence by a sum of delayed 
and scaled unit-sample sequences, suggests that a linear 
System can be characterized by its unit-sample response. 
Epecitically, det h (n) be the response of the system to 


k 
a@(ni—-k), a unit sample occurring at n=k. Then 


co 


y(n) = TI \ x(k) d(n-k) ] or, 
a, 


k= -00 
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y(n) = x(k) T[d(n-k)] = x(k) ie! 
k= -O0O k= -©O 
Thus the system response can be expressed in terms of the 
response of the system to d(n-k). 

The class of shift invariant systems is characterized 
by the property that if y(n) is the response to x(n), then 
y(n-k) is the response to x(n-k), where k is a positive or 
negative integer. When the index n is associated with time, 
shift-invariance corresponds’ to time-invariance. The 
property of shift invariance implies that if h(n) is the 
response to d(n), then the response to d(n-k) is simply 


h(n-k). Therefore oo 


y(n) = x(k) h(n-k) 
k= -—OO 
and any linear shift-invariant system ics completely 
characterized by its unit-sample response h(n). 
A subclass of linear shift-invariant systems consists 
of those systems for which the input x(n) and the output 
y(n) satisfy an Nth-order linear constant-coefficient 


difference equation of the form 
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If the assumption of causality is made about the 
system, a linear difference equation provides an explicit 
relationship between the input to the system and the output 
of the system. This can be seen by rewriting the previous 


equation as 


N M 
y Ga). i= 3 c y(n-k) + Ds d x(n-r) 
k rr 
k=1 r=0 
where c =-a /a and c =b /a 
k KE 0 Poe ruc O 


Thus the nth value of the output can be computed from 
the nth value of the input and the N and M past values of 
the output and input, respectively. The difference equation 
not only represents the system for theoretical purposes, 
but it may also serve as a computational realization of the 
system. The z-Transform makes use of this property to 


realize systems. 


B. THE z-TRANSFORM 


1. Description of the z-Transform 


The z-transform plays an important role in the 


analysis and representation of discrete-time linear shift- 


invariant systems. The z-transform, X(z), of a sequence 
x(n) is defined as co 
-n 
x(Z) = 5 SO NEA 
n= —-Ww 
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where z is a complex variable. This representation of the 
z-transform is referred to as the two-sided z transform. The 
one sided z-transform consists of the same summation for 
terms of n greater than or equal to zero. For the case 
chat x(n)=0 « forram<O,re the: ones sided)» and twol+ sided 2z 
transforms are equivalent. 

‘By expressing the complex variable z in polar form as 

w 
Z = a » the z-transform can be interpreted as the Fourier 
transform of x(n) multiplied by an exponential sequence. 
Heres = 1, that is’, for {zl = 1, the: z=transform is equal to 
the Fourier tranform of the sequence. 

2. The Bilinear Transformation 
The transfer functions of analog systems are most 


Seften expressed in terms of s'= jw (see section III B.). 


This corresponds to the analog frequency response of the 


system. This analog frequency response can be "mapped", 
that is, tranformed to the z-plane from the s-plane through 
the use of the bilinear transformation. The effect of 


utilizing the bilinear transformation is to convert a system 
transfer function in terms of the variable S into the system 
Pranster function in terms of the variable z, The 


transformation itself is: 


and 
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where T is the sampling period, that is, the time between 
data samples. 

Thus a transform can be made from one plane to the 
other. ingthisgemanner;onther transfer function,” HCz)g00fva 
system may be obtained. 

The bilinear tranformation equations may be shown to 
hold in general, and the use of this transformation may be 
shown to yield stable digital filters from stable analog 
felters [Ref. 7]. The bilinear transformation maps the 
imaginary jw axis in the s-plane onto a unit circle (of the 
region of convergence) in the z-plane, with the left half 
s-plane mapped onto the region inside the circle and _ the 
right hand (region of instability) s-plane mapped onto the 
region outside this circle [Ref. 8]. A complete discussion 
of the z-transform is available in several texts, some of 


which are listed in the Bibliography. 


C. THE DIGITAL COMPUTATIONAL ALGORITHM 

In implementing a digital filter on a digital computer 
such as the IBM 3033, the input-output relationship of the 
signals through the system being synthesized must be 
converted to a computational algorithn. The algorithm is 
Specitied in terms of aset of basic computations of 


elements. For the implementation of discrete-time systems 
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described by linear constant coefficient difference 
equations, such as the AN/ASQ-81, it is convenient to choose 
as these elements the basic operations of addition, delay, 
and multiplication by a constant. The computational 
algorithm for implementing the filter is then defined by a 


structure or network consisting of an interconnection of 


these basic operations. For a system transfer function of 
the form M 
=.2 
Deez: 
k 
k=0 ¥(z) 
H(z) SS see ee eee-=-— = —=— = 
N ECZ) 
l- Ne az 
k 
k=] 


the difference equation relating input and output is easily 
written down directly from the system function and is given 


by M 


N 
y(n) = 3 a y(n - k) + ‘3 b x(n - k) Pref. 7] 
k k 


=1 k=0 


oe 


This difference equation can be interpreted directly as a 
computational algorithm in which the delayed values of the 
input are multiplied by the coefficients b , the delayed 
values of the output are multiplied by me coefficients 
a and the resulting products are added. It is now easy 


to see the process to be followed in obtaining the 


computational algorithm for the AN/ASQ-81 magnetometer 
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transfer function. The z-transform of the system transfer 
function is obtained through the use of the bilinear 
transformation, and is then converted into a difference 
equation relating input and output signals, thence to a 
FORTRAN computer program. A table of z-transforms of system 
functions is included in Appendix B. 

In the FORTRAN computer program realization of the total 
system computational algorithm, each filter ‘block is 
transformed into a separate difference equation and 
algorithm. This was done to enable a "building block" type 
approach to the program, and to minimize computational and 


roundoff errors. 


D. THE CASCADE FORM OF THE COMPUTATIONAL ALGORITHM 

Even though the direct form realization of the digital 
filter design may be perfectly satisfactory in a theoretical 
sense, it) may be fess than desirablekin the’scontext of 
realization through the use of a general purpose computer of 
fixed register length. The parameters of a digital filter 
are usually obtained with a high degree of accuracy, which 
Results in a faithful realization of the desired systen. 
When these parameters are quantized, as in a finite memory 
register within a computer, the frequency response of the 
resulting digital filter may differ appreciably from the 


original design. in fact, the quantized filter may fail to 
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meet design specifications although the unquantized filter 
does. [Ref. 7] 

The sensitivity of the filter response to errors in the 
filter parameters is dependent upon the structure of the 
fileer realization. Therefore, in the event of an 
unacceptable change in the frequency response of the filter 
due to quantization errors, it is often possible to minimize 
the effect of these errors through an alternate filter 
realization structure. An alternate structure to the 
previoulsy discussed direct form realization is the cascade 
form realization. 

The direct form network structures were obtained directly 
prom the system.function H(z). written in the form of a ratio 
of sums. tf this) ratio is) factored anto a product of 
polynomials of the form 


[(N4+1)/2] i = 


H(z) = A | |  e---------------------- 


a Z 
ee 1k 2k 


this product represents a general distribution of poles and 
zeros and suggests a set of structures consisting of a 
cascade of first and second-order subsystems. There is 
considerable freedom in the design of the subsystems, but it 
is best to realize the systems using a minimum of storage. 
ine expression of H(z) in this form indicates the 


presence of poles and zeros in pairs. If poles and zeros 
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are not present in pairs, one of the coefficients B or a 
2k 2k 
will be zero as appropriate. An implementation of such a 
cascade structure with the use of minimum memory can _ be 
obtained through a direct form II realization of each second 
order subsystem using techniques similar to the direct form 
implementation utilized previously. A cascade realization 


of a sixth-order system, such as the ASQ-81 system, using a 


direct form II realization of each second order subsystem 


would appear as in Figure 4.2 below. 


Paigure 4.23 Cascade Structure With a Direct Form II 
Realization of Each Second Order Subsystem, 
[Refs 7] 


There is, theoretically, considerable flexibility in the 
Manner in which the poles and zeros are paired together and 
in the order in which the resulting second-order subsystems 
are cascaded. However, although all such pairings and 
orderings are equivalent for infinite-precision arithmetic, 
they may differ considerably in practice owing to finite 


word length effects of roundoff and truncation. 
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Pagure 5:1; Output of First Stage Filter: of Digital Filter 
Computer Program With Sinusoidal Input in Simulation. 
Uatortunately;* the second® stage outputeof the filter 

showed an instability within the program design, indicated 

by the output of the filter being a sinusoid of increasing 


Magnitude, as indicated in Figure 5.2 below. 


(units) 


Amplitude 
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Ereure 5.2: Output of Second Stage Filter Design With Input 
orial Sinuseuds 
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The stability of the third stage of the filter design 
was investigated by inputting the sinusoid directly to the 
maerd yfilter,,. and,» .found.to’be stable. A @beck,"of «the 
derivation of the equations, coefficients, and programming 
steps of the second (unstable) filter of the design failed 
to indicate the cause of the instability. 

Computation of the poles of the z transfer function, 
H(z), of the second stage of the filter confirmed the 
instability of the design. The poles were computed to be: 
meet 8.82161, | 1.07 £ 0.13401, | O.SG11,, and’ 1.1564. Of 
these six poles, three lie outside the region of convergence 
for the z-plane, that is, within the unit circle discussed 
previously in Chapter IV. 

The second) stage! of +) the, filter: was therefore 
redesigned using the cascade form of the direct form 
Bealization (direct form II), and tested in simulation. A 
copy of the software used in the simulation is enclosed in 
Appendix F. 

The output of all three filter stages of the program 
were stable, as indicated in Figures 5.3 through 5.7 below. 
The amplitude decrease and phase shift expected were 
observed. The "damped overshoot" of the second stage output 
Hm aue to the fact that, for values of the input function 
prior to time zero in the simulation, utilized in the input- 
output signal difference equations for the filter, the input 


Signal was. set at 0. This resulted in an instantaneous 
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Exange of the “input signal from 0 “to the finite value 
introduced in the simulation at time 0+. The "overshoot" of 
mre. filter is the (filter's attempt to "match" this 
instantaneous jump in magnitude of the input signal. When 
the input signal to the filter in the simulation is zero 


at time zero, this overshoot effect does not occur. 


TE 
a 
HL 


005s 


000 


(tenths of units) 


-005 


Amplitude 


-010 


Papure 5.3: Input) Signal to Digital Filter Program. A 
Sinusoid of Frequency 0.1 HZ and Amplitude + 1. 
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Pigure 5.6: Output of Third Stage of Digital Filter Program. 


This simulation was run with inputs of sinusoids of 
Warzo0us frequencies in order to check the stability of the 
filter design at frequencies throughout the operating range 
of the AN/ASQ-81 magnetometer. In all cases, the design was 
stable, and the expected amplitude changes and phase shifts 
occured. 

2. Noiselike Inputs 

The simulation was also run with inputs consisting of 
a sinusoid of a frequency which should be passed through the 
AN/ASQ-81 added to sinusoids of frequencies which should 
have been filtered by the magnetometer and random noise. The 
filter performed as expected, with the sinusoid of a 
passable frequency passed by the filter, and spurious noise 
and sinusoids attenuated severely. The sresubtes ‘of a 


simulation consisting of a sinusoid of passable frequency, a 
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filterable sinusoid, and uniformly distributed random noise, 
all of amplitude + 1, are presented in Figures 5.7 through 


5.10 below. 
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figume S,o. Output of Farst Filter” Stage’. 


oy 


— 


hues vase Seaee 


7 


> 


_—— 
—_—a. 
acento: § 
—_ ee. 
= ——- 2D be ‘ 
7s 
_ ~ 
a 


ae Oe - “ 


00e 


Amplitude (units) 


-002 


ee ee ee a 


(units) 


000 


Amplitude 


-001 


—_—— a aa a a a a ae a a a ee ae ae ae ia es i ee ee ee ee ee 


Figure S710: Final Filter) Stage Output. 


As can be seen, the digital filter program succeeds 
in filtering out random noise and signals of frequency 
components above the band pass of the magnetometer. 

In order to ensure that the digital filter 


representation of the magnetometer has the same amplitude 


38 


a -_————— = = ee ee fe ee 


Pee 22129 Speree SA ito iD “ae F 
7 


vauté' peels Usable ae eae 


A os | 


(estar, Pet 


(3NTUnia) sm ?Y 


2 — i ee ee ee eS ee -——— @ © © Ge @2e@ = eo2 e-= 


“die 4 eb, ag; i) Rowe * uD es 


Ate sehdn actas, sus Perret 


uiJ.30 shag Head ats ayess etoon qe 


son? ay 2% hs ed teh 105 ai 
Jotatengan 2d3- Yo qehzatas 


ae a‘ 


versus frequency characteristics of the AN/ASQ-81 
magnetometer, a simulation program was written which inputs 
sinusoids of varying frequencies and computes the Root Mean 
Square (RMS) value of the filter output and the _ signal 
input, then computes the decibel (dB) attenuation of the 
filter at that frequency. A copy of this program is 
included in Appendix G. A plot was made of the _ dB 


attenuation versus frequency for the filter and compared 
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with the measured frequency performance of the AN/ASQ-81 
magnetometer, which was supplied by Texas Instruments, Inc., 
and does not include the effects of the fixed high pass 
fritter. Consequently, the data shown in Figure 5.11 is a 
comparison of the data supplied by Texas Instruments and the 
output of the test program, which also does not include the 
fixed high pass filter. As can be seen, the performance of 
the filter is extremely similar to that of the magnetometer 


itself. 


3. Anderson Function Simulations 
The next step in the simulation phase was_ the 
introduction to the filter of Anderson function simulations. 
The shape of the signal amplitude of the output of a 
Magnetometer passing through the sphere of influence of a 


Magnetic anomaly (submarine) is a function of the dip angle 


of the geomagnetic field, the magnetic heading of the track 
of the magnetometer (or the aircraft), the Magnetic 
heading of the anomaly (submarine) dipole, and the 


lateral range between the magnetometer (aircraft) and the 
anomaly. Anderson functions [Ref. 9] are mathematical 
representations of three basic components of signals which, 
when taken in various linear combinations, describe the 
shape of these anomaiy signals. The equations for the 


Anderson functions are: 
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(First Anderson Function) if Gla Be) 


WO 
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(velocity)x(time) 


range at CPA 


where 


or, a dimensionless parameter. defined as the distance 
traveled along the magnetometer (aircraft) track divided by 
the slant range at closestpoint of approach (CPA) 
(Second Anderson Function) fy = eB x PE 
(Third Anderson Function) He = Rx e = ac x ne 
2 1 0 

The Anderson functions were introduced into’ the 
filter program in a noise-free signal environment in order 
to observe the output signal and ensure that it was a "MAD- 
like" signal. A rigorous determination of the actual output 


Signal would have been extremely difficult to obtain, so a 


comparison was made with the output of a computer simulation 


program provided 'to NPS) by\ Mr. Joe Rice of Texas 
Instruments. When the sampling rate of the program was 
adjusted to equal that of the Texas Instruments program, 8 
HZ, the two program outputs were observed to be very 
similar . The Anderson function simulation inputs.) and 


outputs of the program are depicted in Figures 5.12 through 
Bas LS. The Texas Instruments program outputs were obtained 
im the form of time series plots of discontinuous data 
points, and were therefore not conducive to replotting for 


comparison. 
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Figure 5.12: First Anderson Function Input. CPA at Time 
2.5 Minutes. 


ee a 


Pigmrere 13% Filter Output for First Anderson Function 
Input. 
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Fagure 5.14: Second Anderson Function Input. CPA at Time 
2.o°Minutes. 


(minutes) 


Figure 5.15: Filter Qutput for Second Anderson Function. 
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Figure 5.16: Third Anderson Function Input. CPA at Time 
259° Minutes. 
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Figure 5.17: Filter Output for Third Anderson Function. 


The filter output forralblwecthree Anderson function 
inputs did appear to be "MAD-like" signals, and did closely 
resemble the Simulation output obtained from Texas 
Instruments, Inc. 

4, Impulse Function Response 
The response of the filter program was also observed 


when the input was a unit impulse function. Again, the 
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output was compared to that of the Texas Instruments' 
computer program. The outputs of the two programs were 
observed to be, again, very similar, as can be seen in 
Figure 5.18, -where the response of the NPGS filter is 
represented by a solid line and that of the Texas 


Instruments filter by a chain-dash line. The abrupt "jumps" 
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of the Texas Instruments response are due to the translation 
of the output plot supplied to this plot. The plot supplied 
by Texas Instruments was, again, discontinuous points o£ 


poor resolution, and it was necessary to interpolate values 
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in order to generate Figure 5.18. This resulted in the 
broken appearance of the plot. Even so, the similarity of 


the outputs can be observed. 


B. EQUIPMENT SETUP 

Following the simulation phase of the experiment, actual 
magnetic field measurements were introduced to the filter in 
order to test the response of the filter. Magnetic field 
measurements were made at the La Mesa field test site near 
the Naval Postgraduate School in Monterey. The output of an 
AN/ASQ-81 magnetometer, a Schonstedt magnetic field sensor, 
and the school's coil sensor, oriented along the Earth's 
magnetic field, were pulse code modulated (PCM) and 
transmitted via VHF radio to recording devices at_ the 
Postgraduate school. The recording of a two hour long data 
collection period was transferred to digital data tape for 
use by the school's IBM3033 general purpose mainframe 
computer, 

fie tne first test of the digital filter program, the 
output of the Schonstedt sensor, which represents 
fluctuations of the Earth's total field, was used as. the 
input to the computer program. A comparison of the output 
of the computer program, with this approximation to the 
total field fluctuations as input, to the output of the 
AN/ASQ-81 should provide an indication of the proper 


functioning of the computer filter program. The results of 
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the test are shown in Figures 5.19 through 5.21 on the pages 
following. Figure 5.19, the Schonstedt sensor output, shows 
several instances of PCM dropouts, that is, occasions where 
the pulse code modulation signal was not correctly read by 
the computer for some reason. At such occurances, the data 
point value used by the computer is a random number and does 
not reflect the true value of the data. The problem with 
these PCM dropouts is that the computer does not recognize 
them as invalid data points and will use them in 
computations. This can (and does) cause problems in the 
computation of Fourier transforms, spectral characteristics, 
etc. Additionally, this will also impact the proper 
funetioning of the digital filter program which is the 
supject of this thesis. PCM dropouts are visible at times 
ooo, 142, and 220 through 226 seconds on the plot of the 
Schonstedt sensor output. An examination of Figure 5.20, 
the filter program output, reveals the programs attempt to 
"follow" these PCM dropouts. It should be recalled that 
previous simulations indicated the filter's tendency to 
Etoiilow sudden changes in the input signal, watil <a 
relaxation time required for the filter to steady out. This 
effect is apparent in the output of the filter program at 
times corresponding to those of the PCM dropouts in the 
Schonstedt sensor's time series plot. It can be seen that 
this overshoot tendency resulted in an output significantly 


different from the actual AN/ASQ-81 output at these times. 
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If the PCM dropout induced differences are neglected, it can 
be seen that the shape of the output of the filter program 
is remarkably similar to that of the AN/ASQ-81 magnetometer, 
although noisier. Note that the output of the AN/ASQ-81 
magnetometer exceded the maximum voltage amplitude which the 
pre-amplifiers of the data collection system were able to 
handle and resulted in a truncated signal from time 40 to 
time 60 seconds. It can still be seen, however, that the 
filter program output is very similar to the time signal 
which would have been displayed without this truncation. 

It should be noted that the amplitudes of the time series 
signals of the program output and the AN/ASQ-81 magnetometer 
differ considerably. In the case of the program output, the 
peak amplitudes are on the order of 1.4 nanoteslas along the 
vertical scale, while the peak amplitudes of the output of 
the AN/ASQ-81 magnetometer are on the order of OL 
Manoteslas along the vertical scale. This is because the 
input signal to the filter program is an approximation to 
the total field difference time series signal, and _ some 
amplitude difference could reasonably be expected. The 
intent of this initial test was to investigate the output 
time series’ shape, and an exact correlation was not 
expected. le is worth noting that the digital filter 
program will perform its function on any time series signal, 


regardless of units. This means that a signal may be 
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operated upon either before or after conversion from 
whatever units it was originally measured to magnetic field 
strength units. 

Therefore it appears that the digital filter program is 
functioning properly. When a close approximation to the 
fluctuations of the total field time series signal is used 
as the input to the computer program, the output of the 
program is similar to the time series output of an AN/ASQ-81 
magnetometer. 

The final stage in the testing process was a conversion 
of the time series output voltage signal of the coil antenna 
sensor, which was aligned along the Earth's magnetic field, 
into a total field fluctuation time series representation 
for the same time period as before, and then to use this as 
Pe input to the digital filter progran. A comparison of 
the resultant time series output of the program with the 
actual AN/ASQ-81 magnetometer output would validate the 
proper functioning of the program. 

Conversion of the time series antenna sensor output 
voltage signal into total field fluctuations in nanoteslas 
was accomplished through the use of a computer program 
designed by Capt. Kurt Stevens, USAF, a student at the Naval 
Postgraduate School, as his Master's thesis [Ref 10]. The 
output voltage time series is stored in an array, then a 
Fourier transform is performed on the stored data, resulting 


in the Fourier spectrum of the data. ihas spectrum is 
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corrected for the characteristics of the coil antenna sensor 
to obtain the Fourier spectrum of the total field data. A 
reverse Fourier transform gives the time series signal for 
total magnetic field in nanoteslas. 

This time series signal was used as the input to. the 
digital filter program and compared with the output of the 
AN/ASQ-81 magnetometer. Figures 5.22 through 5.24 show the 
raw coil antenna data, the total field time series data, and 
the program output time series for a 6 minute period of the 
test. Figure 5.22 shows the raw coil antenna data. series. 
The number of PCM dropouts should be noted, as these will 
influence the performance of the filter program. Figure 
5.23 shows the computed total field time series. Note that 
the PCM dropouts evident on the raw time series plot are 
evident on the computed total field time series plot also, 
and thus inputted to tthe filter program as valid data 
points. Additionally, there are two "jumps" in the plot of 
Botal field  fluctuationmacrigure 5.23) which are also 
inputted to the filter program as valid data points. These 
"jumps" are located at 128 and 256 seconds and are caused by 
the method of processing blocks of data for the conversion 
mo total field fluctuation: A block of 128 seconds of data 
is processed at a time, and the results of each block are 
stored in an array. this results in a siasht amplitude 


difference between the last data point of one block and _ the 
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frst data point of the next block of data. thas) siight 
difference is manifested as a signal jump. 

AS comparison of Figures 5.24 (program output) and 5.21 
(AN/ASQ-81 output) show that the filter program gives a time 
series output very similar to that of the actual 
Magnetometer. The first 20 seconds of the program output is 
somewhat dissimilar to that of the AN/ASQ-81, due either to 
Ptcmoinitial “start up" delay of the filter program or to 
distortion of the total field fluctuation time _ series. 
There is a PCM dropout at time 11 seconds which contributed 
Eos che distortion. 

Following this, however, it can be seen that the program 
output is very similar to that of the magnetometer, except 
at 128 and 256 seconds, which show the effects of the false 
signal jumps caused by the total field fluctuation 
conversion. There is also a noticeable four to five second 
time delay between the AN/ASQ-81 output and that of the 
filter program. As| this delay is’ mot evident in a 
comparison of the AN/ASQ-81 output and that of the filter 
program with the Schonstedt sensor as the input, it can be 
inferred that this time delay is caused either by the 
program which converts the raw coil data to total field 
fluctuation data, or by a phase (and hence time) change of 
pues voltage signal due te the coil sensor itself. A 
comparison of the raw coil data in Figure 5.22 to the 


converted coil data in Figure 5.23 indicates no time shift, 
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and hence the deduction can be made that there is a time 
delay inherent within the coil sensor itself. 

Other than the differences of the four to five second 
time delay and the distortions caused by the false signal 
jumps, the program output is extremely similar to the output 


of the AN/ASQ-81 magnetometer. 
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VI. CONCLUSIONS 

The intent of this thesis was to design and test a 
digital filter computer program which would, when given a 
time series input of fluctuations in the total magnetic 
field, deliver an output time series representation of the 
output of an AN/ASQ-81 magnetometer. This purpose has been 
realized. 

The computer program contained in Appendix I has_ been 
proven to output a time series signal which is very similar 
to that of the magnetometer. The major limitations of the 
output signal are a finite time delay of about five seconds 
between the AN/ASQ-81 magnetometer signal and the output 
signal of the program, a sensitivity of the program to false 
data points such as those caused by PCM dropouts and false 
signal jumps caused by processing large data blocks, and the 
inherent limitations of the program caused by its dependence 
on the use of digital data tapes and the IBM 3033 mainframe 
computer. 

The five second time delay is not considered to be an 
important limitation to the program, as it was intended as a 
research tool for programs currently in progress at_ the 
Naval Postgraduate School. Instances where this time delay 
might become important would be in areas of simultaneous 
comparison of the program output signal with an actual 


magnetometer, in target location algorithms using time 
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delays, or in correlation studies between different sensors. 
In correlation studies using coil sensors, the effects of 
the time delay would cancel out, as all coil outputs would 
be similarly delayed. In target location algorithms, the 
target location errors due to the time delay could be 
adjusted for simply while in computer simulation, and flight 
testing could not reasonably be accomplished without the use 
of an actual magnetometer as the sensor. Lastiy, in .a 
comparison of the program output with an actual sensor, the 
time delay can, again, be compensated for. In short, these 
limitations are not considered excessive, especially as the 
apparent cause for the delay is not the filter program. 

In the primarily intended purpose of the filter progran, 
Magnetic noise studies, the time delay is not considered to 
be a problen. 

The problem of false data points caused by PCM dropouts 
and signal jumps due to conversion to total Eield 
fluctuations is more_ serious. False data points cause 
inaccuracies in the output time series and could adversely 
aerect later projects. Unfortunately the PCM dropout 
problem is one which is endemic to the data collection 
system presently being used at the postgraduate school, and 
not to the filter program itself. It is imperative that 
users of this program are aware of the PCM dropout problem 


ene Of the effects it may entail upon their specific 
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research, A large number of PCM dropouts in a time. series 
could render that series unusable. Similarly the case of 
the false data jumps caused by conversion to total field 
figectuations is not within the filter program. Further 
investigation of this problem is necessary in order to 
eliminate it. 

The last problem, that of reliance upon the digital data 
tape/IBM 3033 computer system, is, like the PCM problem, one 
which is not endimic to the filter program but rather to the 
data collection system being used. A change of data 
collection system may, at some future time, remove’ the 
reliance upon the PCM/digital tape/IBM 3033 data system (and 
hence too the data block conversion problem which results in 
falso data jumps), but this is unlikely at this time. Users 
should be aware of this dependence and of possible effects 


upon specific research projects. 
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APPENDIX A 


AN/ASQ-81 FILTER TRANSFER FUNCTIONS 


Fixed High Pass Transfer Function; 


Selectable High Pass Transfer Functions: 


A. 0.04 HZ H(S) = 
2 


40.82834 §S 45.283170S 


SOS I 2) ee ES ork aa >. 4 Pee 


0762834 S$. + 12.52096 § + 1 aD. ZOo17 (S 
2 
57.576688 § 


ee a a ea ea ee = 


Pdi 00999 3S + 1 


ac ae ee a eee a (ee fe ees ee ee ie 
BPS a6688-Si + 7 h1498.8 4 1 
mee 0.06 HZ. H(S) = 
2 2 
18.2459] s 20.12587 § 


25,96904 S + 4.94332 S + 1 


@ 0.08 HZ H(S) = 
2 


2 


10.20708 §S 1153208075 


WEBOGl7 S- + 3.70749 S + 1 
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Selectable Low Pass Transfer Functions 
eee). 2 HZ H(S) = 

1. 1 
ee cee ee ee ee ee es we we ee ee ee ee ee ee we ee we we we we we we we ee y'< Mian MOR paar 
055143 s 2 Phi sOV4GL eS «tal 0.2501 S- = 0.6209 S + 1 
Be 0.4 HZ BCs) = 

1 1 
_ 2a a ea a ee! a ee) Cee ae ae a xX 2 Rs Tae ©, PS ee eee eee 
Seo7e56 S~ + 0.53706 S + 1 O306252°S> + 0531044 S$ + 1 
re 0). 6 «HZ H(S) = 

1 1 
a= >> >> >>> X  --------- aaa ae a 
0.03492 S' + 0.35804 S§ +1 0.02779 sS. +20920696 S + 1 
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APPENDIX B 


AN/ASQ-81 Z TRANSFORM FILTER TRANSFER FUNCTIONS FOR DIRECT 
FORM I REALIZATION 


Por fixed high*pass filter: 


f"="APHPISZ." = AFAR IZ" 


where BFHPO, BFHP1, BFHP2, AFHP1, . AFHP2 are constants 

tabulated in Appendix D. 

For selectable high pass filter: 
-1 


a) =3 
BSHPO + BSHP1*Z + BSHP2*Z + BSHP3*Z 


ih 


+BSHP4*Z~* + BSHP5*Z~- + BSHP6*z ~° 
H(Z) = noo 7a tac a ae UR gw poy Fae a De Ae ee SO a Too ao 
1 - ASHP1*Z  - ASHP2*Z ~ - ASHP3*Z 
- ASHP4*Z~* — ASHP5*Z~> - ASHP6*z~° 


where, for low frequency cutoff of 0.04 HZ: 


BSHPO = 0.99471378 

BSHP1 = -5.9682827 ASHP1 = 5.9894021 
BsurP2 = 14.920707 ASHP2 = -14.947051 
BSHP3 = -19.894276 ASHP3 = 192894225 
BsHrp4 = 14:;920707 ASHP4 = -14.894327 
BSHP5 = -5.9682817 ASHP5 = 52907 2141 
BSHP6 = Orgsa7 i372 ASHP6 = -0.98945296 
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For low frequency cutoff of 0.06 HZ: 


BSHPO = 059920613 

BSHP1 = -5.9524928 ASHP1 = 5.9841070 
BSHP2 = 14,881232 ASHP2 = -14.920650 
BSHP3 = -19.841643 ASHP3 = 19.841528 
BSHP4 = 14.881232 ASHP4 = -14.,841757 
Beges = =-5.9524927 ASHP5 = 5, 9209919 
BSHP6 = 0-992038212 ASHP6 = -0.98422128 
mom tow frequency cutoff of 0.08 HZ: 

BSHPO = 0.98945806 

BSHP1 = -5.9367483 ASHP1 = 5.9788144 
BSHP2 = 14.841871 ASHP2 = -14.894276 
BSHP3 = -19.789161 ASHP3 = J19..788959 
BSHP4 = 14.841871 ASHP4 = -14.789365 
BSHP5 = -5.9367476 ASHP5 = 5.8948841 
BSHP6 = 0.98945802 ASHP6 = -0.97901720 
wom low frequency cutoff of 0.1 HZ: 

BSHPO = 0.98684156 

BSHP1 = -5.9210493 ASHP1 = 5.9735244 
BoHP2 = 14.802623 ASHP2 = -14.867941 
BSHP3 = -19.736831 ASHP3 = 192730516 
BSHP4 = 14.802623 ASHP4 = -14.737149 
BSHP5 = -5.9210491 ASHP5 = 5. 8688889 
BSHP6 = 0.9868415 ASHP6 = -0.97384065 
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For selectable low pass filter: 


BSLPO + BSLP1*Z + BSLP2*Z + BSLP3*Z + BSLP4*Z 
J = ASLP1*Z— = ASLP2*7 > = ASLP3¥z7 -° — aSLP4*z ~* 
mhere BSLPO, BSLP1,. BSLP2, -BSLP3, -BSLP4, ASLP1, ASLP2, 


ASLP3, ASLP4 are constants tabulated in Appendix D. 
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APPENDIX C 


AN/ASQ-81 Z TRANSFORM FILTER TRANSFER FUNCTIONS DIRECT 
FORM II REALIZATION 


Por fixed high pass filter: 


1 - AFHP1*Z~+ - aFHP2*Z~* 
where BFHPO, BFHP1, BFHP2, AFHP1, AFHP2 are constants 


tabulated in Appendix D. 


or selectable high pass filter: 


= =2 
be-o2¥Zut +fZ 
H(Z) = ASHP] x -------------5-----------5 x 
Jes ASHPS*Zde1<; ASHP2*Z 
=i £2 24 22 
me 2*2 + 2 W 12470 VENZ 
eo sas = == SS ssen aeas x Satan UNE REE aoe Te > 
fe] ASHP4*Z.~ — ASHP5*Z 1 = ASHP6*Z ~~ = ASHP7*Z 


mere ASHP!, ASHP2Z, ASHP3, ASHP4, ASHP5, ASHP6, ASHP? are 
constants and tabulated in Appendix D. 
For selectable low pass filter: 


Pal a2 43 yh 
BSLPO + BSLP1*Z + BSLP2*Z + BSLP3*Z + BSLP4*Z 


i ] ASEPI*Z - = AGLP2™Z --= ASLP3*Z -— ASLPA*Z ~~ 


waere) BOLPO,.,’ BSLP1,'-BSLP2;° BSLP3, BSLP4, ASLP1, ASLP2, 


ASLP3, ASLP4 are constants tabulated in Appendix D. 
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APPENDIX D 


Z TRANSFORM REALIZATION DIFFERENCE EQUATIONS 


With reference to Figures 3.2 and 4.2, the following 
difference equations are used to model the AN/ASQ-81 
Magnetometer filter transfer functions. The input to the 
eexeag high pass filter is called SIG(I), where I is the 
current data sample. The output of the fixed high pass 
filter, which is the input to the selectable high pass 


Bolter, is YO(1), and the output of the selectable high 


pass filter, the input to the selectable low pass filter, 
ieeecatted YPO(I). The output of the filter is called 
ASOCT). (I-1) denotes a time delay of one sample, and so 


forth, and the symbol * denotes multiplication. 
Mer the fixed high pass filter: 


YO(I)=BFHPO*SIG(I) + BFHP1*SIG(I-1) + BFHP2*SIG(I-2) + 
AFHP1*YO(I-1) + AFHP2*YO(I-2) 


where: 

BFHPO = 0.9980499222938581 

BFHP1 = -1.9960998445877161 AFHP1 = 1.9960983216843922 
BFHP2 = 0.998049922238581 AFHP2 = -0.9961013674910398 


For the selectable high pass filters: 
XI(I) = ASHP1*YO(I) + ASHP2*XI(I-2) + ASHP3*XI(1I-3) 
KXII(I) = XI(I) + XI(I-2) - 2*XI(I-1) 
metre) = XTI(L) + ASHP4*XTII( 1-1) + ASHPS*XIII(1I-2) 


MEVCT) = XITIC(L) = 2*XTETCI-1) + XILICI-2) 
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a ae 


XV(I) = 


TPOCT) 


For the 


ASHP1 = 


ASHP3 


ASHP5 
ASHP7 = 
For the 
ASHP1 = 


ASHP3 


ASHP5 = 
ASHP7 = 
For the 
ASHP1 = 


ASHP 3 


ASHP5 = 
ASHP7 = 
For the 
ASHP1 = 


ASHP3 


ASHP5 


ASHP7 = 


For the 


ASQ(T) 


XIV(I) + ASHP6*XV(I-1) + ASHP7*XV(I-2) 
= XV(I) - 2*XV(I-1) + XV(I-2) 
low frequency cutoff at 0.04 HZ: 


0.994713789347288 ASHP2 -0.9952196910157882 


1299 352157250322473 ASHP4 1..9962028201103847 


-0.9962082013015601 ASHP6 1.9979855321466768 
-0.9979897681491607 

low frequency cutoff at 0.06 HZ: 

0. 9920821277199399 ASHP2 = -0.9928381306174365 


1.9928247245317958 ASHP4 = 1.9943056083414792 


-0.9943177045263504 ASHP6 = 1.9969766455640039 
-0.9969861717656565 
low frequency cutoff at 0.08 HZ: 


0.9894580558875787 ASHP2 -0.9904622611337722 


1.9904384565912725 ASHP4 1 .9924092959984783 


-0.9924307799269113 ASHP6 1.9959666590811389 
-0.9959835860201267 
low frequency cutoffl at 0.10 HZ: 


0.9869415560096681 ASHP2 -0.9880929619779491 


1.9880549117849344 ASHP4 = 1.9905139074397893 


-0.9905474442556225 ASHP6 1.9949555824611562 
-0.9949820174650785 
selectable low pass filters: 


ASLP1*ASQ(I-1) + ASLP2*ASQ(I-2) + ASLP3*ASQ(I-3) 


+ ASLP4*ASQ(I-4) + BSLPO*YPO(I) + BSLP1*YPO(I-1) 


+ BSLP2*YPO(I-2) + BSLP3*YPO(I-3) + BSLP4*YPO(I-4) 
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high trequency. cutoff .at..0...2 


0.0000000452616229 
0.0000001810464917 
0.0000002715697375 
0.0000001810464917 


0.0000000452616229 


high frequency cutoff at 0.4 


.0000006918001209 
.0000027672004837 
-0000041508007256 
.0000027672004837 


.0000006918001209 


high frequency cutoff at 0.6 


0.0000033463317975 
0.0000133853271900 
0.0000200779907850 
0.0000133853271900 
0.0000033463317975 


ASLP1 
ASLP2 
ASLP3 


ASLP4 


ASLP1 
ASLP2 
ASLP3 


ASLP4 


ASLP1 
ASLP2 
ASLP3 


ASLP4 
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3.9082436339591027 
-5,7285022156249328 
3. 7321935213310065 
-0.9119356638511430 
HZ: 
3.8173771378993420 
-5.4670046844062743 
3.4812554457127576 
-0.8316389680077603 


HZ: 


3.7274299052305002 
=5.2152772583906819 
3.2462216641520829 


-0.7584278523006615 
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APPENDIX E 


DIGITAL SOFTWARE FOR SIMULATION - DIRECT FORM (SINUSOIDS AS INPUT) 
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APPENDIX F 
DIGITAL SOFTWARE FOR SIMULATION - CASCADE FORM (SINUSOIDS AS INPUT) 


ooqo0o0o0 000 3000000000090 00 90 90090 0090 09D 000009099" °9°99D 0090909 
DDORSN (APA DOP DDO ANNO FILN DO Re DDO SANS AOE DRDOANAPLIVOM DRDORAIOSLN 
DDD hE hb ON NI OI NII IN GN NY OV. FY CER COD. COVEY COV. OF) So SF Se SP SS SSP SSP LALLA LAL 
AN NAINNAINAINOAIOAIN NAIA AYN OANONIN NAN NINN NON ONIN NNO ON ONIOINOINN OY 
AdSDODDVDOQSCDOV OBS SNS T.o 99 0N 90g aD Va NN 9D3S9N900090NONST9NTDB00009NO 
QAaagaadaagaanAaaaaaaaaaaaaaaaaAaaadaaaaanaadaaaaaaaada 
QWAAAAAAAAAAAAaAQWAGQAAAAAaAA OOOO OOo aoa aaa ea 
Qe tot eet a Sat es I os ef ent gf fd I est’ QI tara 


~~? 
fa qe Be) uJ o - ~ r . 
LreQ Dr ~ - s = ~@ 2 
ZIt~> — ~ 
~JjWOWY — —= 
ap A Pra Sw S) - 
Re=rDw — cen) 2) “4 a; Qs G7 Gs a: 
— ee) a4 S SE fp -e ee 
LIN Sus fan) Bee of Mm Wm OY Ww 
we “uw aw ~ e a rr bt +S + 
eye) (es (a Pa > + tt 
rq VEU at x< Oo wn 2 - -= 
Ue Le Lem) fag = ey tL oO - -e @ = 
5 | u uy Go = e - = os = - - 
dur <I>» oOo Oo WwW = 8 - 
QO moO © > ) = (ee +r a 
Lo> ma mn oa a ° 
Ce Mud ~ oO - > = — j=) 
EF vOO = wm~w Quan Ss — N 
22r0™w ie.) = > - ul Se aed Ee aM 
ORK NnNO @—lwW “~~ — Oi Pade (tol = Ga ed -—~, 
_ od ~w — x«< eye nN FO =m ZF OO YB eN 
eo TD er {Te = ~- qe Fr Oe 2 - 2 = On 
OnaxYO (9) -_=- — 70 Aa Qe -_ © - zm Ne 
<< =) ty 6 <tuy So Nz4O =e -= ese = Ww Oe 
300M Ie 9 Daen ND OH YN DW ia 
WoO mae aed Oo 90O fen =< Ow ~~ nD Wn a Ne 
Lireteo 4 On LI ee & me NY ”n ~ #0 
KA 2 az ~~ OHO LEONRF D - ~ Ue eS 
A 2 za Da —= este FOOD Aw DS 2 hj 
ReetSe=z OF Cm ses AAA &Fe DOD FO we) YA +++ 
A wtoOo ~-eK j= — AQ eD LV 2 D> ez a. —_ ”n es 
WWke—O Wee - <~ TeDd NNMOwH O KF | = (2) zt oe 
on} ZY Caz = - = “Lore tthe . = = o Lan} a n~ 
" tw>eEzt iw Ome Dat ee nme = = & o@ UL, —~ 
Y eAe) Cp) os ase) S00 est MmMNe © © = @ WE x= +e 
WY Ree ~—at ooo oMOt VOU uw (a) ee 
a q et OwRFEr ak MoO Aaefe aaa S — AHO 
— OOwWLOrF NO ~nmy Len Ciy ws) = say 
YQ Ui ee ee Owe ULL NNW WwW Ww Ww ~O 
ry Fas © = 2) ONOm Dee AIOrF KF§ KH KF &- YY aS Qo =™/_ 
- NOCO = - tae ~— LU eye WwW Ww Ww WwW We Ww oN 
Zot = wl >< NW NNAS 5D 3D 5D 2D ww <O <x Nr 
= QnOZOY OFT =ee Aee,- FOOL =. © ft rn" (8) — 1% 
~ TSWOReAyU =e Ome THA wade 2 = 2 = NLNO WwilW ue a 
= Ace we eXolo) (jal=\— haa aS ye bess) ey) maa) O+ 
Pe) any Yn O00 te enunnn- = = = = N22] zw Oe 
jo) YN Ie — maooO ROSIN NN NN ON ao ac QO =a 
a OC <T LU e ~ONM sOOm ee est SS Ht 4 Vee SU ue. ~n 
ml wae VY) = (Swe A @ OY SAH O~wN ewww ™ tO OO Ne 
= ZrPUOwUs swe SOnOLAMY -ovute De Ye Qe We @ re nN w| 
Sat <ae ZF ..£<D A>-<XWLOaOMnaanas 4 as a 4 Neto n  % e@ 
Wa Candee Wteestiyje —-H —- KH KH Sew Aw zee 
FOO OL 2 ak ZSP2ZN AIDING ZL 1 we 
sty FO Oe, Qunuo0erst929N0SI Or KF KF KF KF QONA aa oo 
~—O CI >O = SRY <n Pte) 1 | oat] 
x2 Atkee AaK<t NNVWODDODDDDDDMD D0 DOD DOD D wow Wt wm I 
oK— KWo DDC SPANZAHHHwHeweHe eee H HF iH a 2WYW nN Ww eeu 
OaWN NeZz>Y CS a es ee ed Ss od ae ew aw 
TU met iene ROO OES EAI I tt fcte Pe dae Se te IRD LZ aw a TT 
So Tommi LLU LU LU LULL Le Le ee KOs dt iO 
woe RFOZOK NOHWQ OOO CLCLCLKCLCLYLODLMOLDLDLOLCOQ OF KF YO sta 
RW aA A A A BF 
LI >< ae 
SwWO 
as UE 
~ne™ 


NS NOOOULDOULUUUO Ou YOU 


76 


(TO Bs ETOPURE) ae 


y= * ate 
cVLTeoesacc 
ee eet, tll er 


Da ees a | ie 
oceans oi 
PPT eS oe aA son vee 
pap P= + o = 
ay 7 . rm 


_ 


wh 


a_¥Y¥ 


Eee 
= 4 


= 
Cc 
© 
ec 
fe 
~~ 
= 
= 
° 
ad -_- 
© 
® ai @ a] ’ y 
’ 7 
- _ i © Al —= _ 
7 ¥ id r- | 
~ ~ rr 2 @ @ -— Ge 
= Tt oy 7 
- ke . rm Ww *§ £© fF c. 
> » us = - ie zl = = ~ 7 
* = i-ie : orn 
- e a «= Pe > 2 . C 
: cope: Se of ee, ww 
2 * e =, ¥ ‘ a 
-~ av: gy 
ys its fs ." 
- Si- * reese S&S — 
‘ - “ u - 24 a a = = 
«a ners ir & - > 
o = a > eo e | > 
ont cn) us it fe ond y hed ¢€) 
. = nd —— 5 _— —_ 
e iw ee * = - wrtrs 
- ° r “ — a _— ? 
& om ( 
s . x - pe Saal — 
iy > > Bae 
ma) - e>€ s 
oe he - A é 
We taw is oa y ° ——€ 
“eo cy fy 
= v 
5 u 
; . é «a 7 
o hoe - » 
a & @ ' ‘ 
. = ‘. rT. 
<« _— 
he 9 . ee 
. =. pa 
‘ Sei) Sivas + 
nm i rt ‘7 
cut ae = Ts. i 
_ ~, > e -_ f fe 
=r > Po 
Be Gh fo Cet Orie we 
. Nl rete Bee 
- ipa La —_— width): Baa 
a = . wi we ahead 
a. a A b wal alls en 2 
GD sah Dis a ee Cc 2m 


Can diab i am 


n 


QODVSSSOSQS SSN OD OOSNNS sO BODO 9NIBV 990 N0DNBV 90 VDONV90NSC8NSsIN909N9N00 
DEH DDD AMOS LIVE DDIANOFAN OM DDOQANM- FN OM DHOANAPADY DHOANY 
WL LVN. 00.9 0000 ORE EEE EYRE El DDD DNDDDDDADAARAAIAAVWARAHOSVOS 
NAUNNNNINNANNNNAN NAN NNNANNANN NO NOINIOAINNNNANN ANN COD COOH 
SOV CcOOSCOCOOOVOVNBVIIV NIN DION DNV NSO VASA See 920920000 09D90909N90930 
AGQaAAAAGAAAAAAAAAAAAAAAAAAAAAMAAAAaAAaAAAaAAAAAAAAAA 
2aaaaaaaaaaadaadiadaaadaaaaagaaadgaaaagaadaagaaaaaadadaa, 
ob OO OO Oo oO oe Oe ee 


=~ yw 
~ ON 
enie 
= <4 
~\|~=— 
Wout 
onan 
aAaa, 
2EI525 
Lue i 
fo ayeeyea) 


GE HIGH PASS FILTER 


OUOU 


® ee ee e 

ttt +r + 

STS OS 

-_ =e —_=—_ -_= 

N ANN NN SN 

<P Ts + we wth + 
OT + DO Hib e ath 3 Ot 
aM O & RK NREOIe 
N DoD QO ~AD~w~w~ wewN~ 
2 N Ni 4 ee me 
ee YW OQNDANOLAIWL 
O+iAn 2D + ie e+e Ht + 
TOT sm O of 8&8 a 2 ow @ 


SNOAMNOODNR UNE NNE OI 
DO AIO GON Re ee er wer SS 
PDO NAIDIN & 3h bh hn Ht 
SO oD oD 0% OH % OF HLH 
NO OW? St et POW + UW 
UPFOPHUVG of + of # of 
em ONSTR ON 0 oN © eC! 0 
N ers @ 0 Qe fl ania Ss 
Settee nn 
WU WW at et ae 
TovdwWeTAOVOWu StI 


SELECTABLE HIGH 


E 


ENT FOR TH 


WASHP41" 
ILTER WITH 


S 


mu 


OUOU 


ria) 


~SOOWIL Im 
0 we we ee ees es ee 
en JL Ome Hk fe Se a 
it | | | 
AINO Pu Of 
RS ike en eS oe ee oe 
aaaaaad 
eke aloats leek eis ees 
NNNNNNY 
tactic 


e 
wr 
~~ 
NE eND em 
DOW OPN 
tN tN 
MTFONONY 
CGoo0norade 
oD 2 OD owt 
OWNGO0O# # AHH 
NONNAN ST ect eo 
eeeeoe®? efron 
soon Hou ow it 
TN | ee 
aAnoowLtooudu 
IOOOWWTOOOW 


ae FOR SELECTABLE LOW PASS FILTER WiTH UPPER FREQ. 
03492 
O2779 
(T** 2)/4. 
DDI 


OF O. 


YOUU 


.#] 
= 


er és a or ) 35 
ev ats pa ge 30 it 1en.% i 
Hob 240NTI9237 oof AA Wei?) pyvau 54” 


fe 


“lf 24 ICH 


Ve 


QBVoo0DC ODDO IVI POO VO NON SOV OV COO OSO VSO SN DOO ONV OOO ONOOCV900Sg90 
FDO DRO MSV S AOE DHRAOAAUAITOAOHPODADSANMPNBGH DA DANOSON OF OHO 
DOD DDD. ah BAIN IIA OI NI AION FF OF) COV EFF CFV. C9. SP SP St se SS SSP SLY 
VDDD ADV AAO DA AH HAMA AMD AAA AHM HAMHAMHAOHMAHAMAMHDIOAMAMAD ANA QAMNHcodH 
o0090g90 000000 09000090 00 09 000900090090 000090 00900 3090 090900900 
Qaadaadmadtaaraaaaaadadaaaaaaaaadaadaaaaanadaaaagaqgaaaa 
BAaAAaAaAAARAAARAAAGAAAAAAAAAAAAAAAAAAAAAGAAAAaAAAAAAAAAA 
tot VS Pe PII es IIe i i et eq IQ et ttarr—zt 


INITIALIZE STORAGE REGISTERS 


EQUATIONS FOR SIMULATED SIGNAL 


— 
= —_ 
> i) 
a ] ane 
+ aur | 
fn] + 
(as) Ce | 
> ) oO 
es ras } 
aN _—— = 
“-a— =“ - 
2-4-7 ar) a) = 
+ 3+ ir Lon bee ) Lad Eel 
oe ae Lae % Oe za = 
7. £8 COmO OD — ~ 
~ SNR NNR SS emwOw = oO — 
= He NON ~ 
NNN Nada te es) Ty i 
+ He POM YH ORAS OL ad ad = 
emt Hm tI MAHOU OL +e 3+ (om) 
kr ott omit eH OBOWUeS OO — ee J 
we NN wwe er st IHN GE HHO OU uJ NN Uu 
HN He RH Tee ORM VY % + + 
OPmOW PL TEN ew OMOO~ 19 i ~— 
ONOWNL ++ een DOON 2 “ AQ ao 
ee ee ee ee es ee) fa UL te Ww 
tH © ot eC Joes eH OO uw WLS ttt et tO LU 
NENNE NYS OWuLLh 5D CO —_ NNO Vou hon 
Nw NN wrt te POUL Ze) eAAao ANMTND eO 
eet eel att, 2 le Le Le Lee ee ee: ie OO = O~e~w) 1 ht eg dg ee 
OOOWLiws Soa ogtIMO oO -4 Us PULL Set ret tet IL 
CO de DUS +} LW ww ee <I IMU woe ee? et sno ee ee? osm 
ke*el+elititi~w~wwww C2 DO 1000 ae ee) pd el el bd Bed ed Ot PO) 
a7 oN 08 oN of} HU et ae nonntod w Wood cad aa soaed nef aed wed 19 
Qe set ARM PORNMS CF Oe t|lmnm—nZ O DHOANMTLVO 2 © © 2 & a= 
OD www we we ww 00000000 YN O92797T7— TF Anwwvorr)~eprrs’_sranatnvodeat 
WT CC CC mm mm Nt et tt et at bag tt Pt er tg 
Ss |e Lm La Lo La Lae ee ee ee ee ee De DwO OOF OC Weare ern eer we SOQ 
LOATH YZUINNNNNNNMNY 4+ OYOANO —K- VNETORNMNTOOLULLLLe 2 
Lota wut ooono BY OR PRPINYD YN OAD HH RRR RSet Rea ft 
fe) 
fe) 
N 
Ouw OuVY Oo 


78 


t? 
a) 
wee 


Vet ye 


Ey ) 
‘shu D> 


= 


TINA SLED 2T RET 


TOREDTECTENL2° ML TW ETE Simewee Mert eyere 


-~ 
-_ ; a 4 
\* ; 
* ~~ 
: e £ 
& = + megane 
- we pe 
. e-e fae 
£ 
ea 
t eel ea 
e @ 7 
ve oe 
i - 
ger hae -= 
yw 4 = e 
Oe. i% ca 
~ a be r a 
COT ot q * i 
ioe tou ge > 
a wooo or ec 
or we fhe ok iy 
~ = yarn 24 Se obmee Oh mem -< 
j 6: 4 arr rer => ft 
—. eer 4 |, ~ 
‘ee ene ES See = 
; | ik — 
7 ae es in femal m2 = poe 
a8 


bie bd - cyest.! 


‘er 


QAODDVDOOVOAASIN SOO ODO SOSNSOV SOON Ng NI ONS SNVOANVNIO OTV9SININSCONOCO9D90390 
NOPD OP DODD ANOS OM DDO ODOM DAODAUOIFIADIMH DHRO.ANMAFOOE ODD 
LAVLAVLALVIA LAL 0.9 000000 OME EEE EEE DDDDDDDDODDAAADAADHAAAN I 
DV DVM DVM OVO DVO] MAVAMD OA OMAMOAOAAAHOIAMDMAM AANA AQAA AMAMAMH AAAI AGAIMAIMAAAY 
Plolelolololelelelelelele)\ele\o)lelplolelelelolelelololelelolaloleleloleloleleloelolelslelelalola) 
AQAAARAAaAAAAAAAAAAAAAAAMAAAAAAAAAAaAAAaAAaAaAAaAAAAAAA, 
AAAAAAAAAAAAAAQAAAAAAaAaAaAAAAAaAAdOAA AAA AAAaAOOAO 
Ail fff fo i I I I eI ot i q@ fi fc ft i te et et 


STGNAL 
005 


THe 


DDS NOISE OF t/- 0625 TO 
F=10 HZ OF MAGe +/- 


EMENT A 
A SINUSOID OF 


STAT 
TO 


= 
+ 
— 
— 
— a | 
= “” 
— od 
~ + 
(=) + 
>» eo) 
+ a 
— ~~ 
a WY 
a a 
up =U + 
YN - 
— + o-— 
o00U.=elC F=hl — 
ZN eS ON = 
DG ee Bo | G~ 
-_ ~~ ~ nO 
N = x = + > 
—~ YN e+ = = NH 
rc r+ MO KEK NN ON 
an FTF +H = ADO 
+ a Q waewawe WA 
— 35) 25 S22 eA) 
mm Lb YW Ams IM 
~ oOo xt Inv +o 
_ + + NEE —+ 
<t- 7 <P NO 
Oo FN RQ HS 
a 1 
we et) me we 
e+ OF Sew kK LHHNOM 
— Kee + TO 
Q DD Hwan is > OO 
He 3 OSS ON He 
eS SP Khe <0 Re 
N Q At X~—wHAwOH 
=~ CFC TF eter> edsay 


az te NN PH XKXNY IO 
ONAOD Ll ree 1 CH 
+—a xo +O 
m—N Ie PAinm +MY 
ON RO) CQ tes OO 
ewe TY wer et + 
OCS Onm-r | N>~—w~—= 
Z2rerwr hPa tT TOM 
tome N= & eb Ht NOHO 
wot 1 ON mene or Pa dw) 
Hw O OM tt te OO + 
eC LOD le KH Se >00}-a 
Or Il Leu I) Ke <0 HO 
Wt —Oowe SH <I 
Wee HO I — == |} MNO 
Dt es me, em et wee et oy <T COO 
mg er LL ed gees at er eer jf WY —J~—e 
OO S welh we Se OO RHING 
Zee OTe Sa aa ON 
aa 
a 


100 CONTINUE 


PLOT GUTPUT 


COMPUTATIONS FINISHED ANU ANSWERS STORED IN ARRAYS, 


no 


ILL BE MAGNITUVE AND 


W 
LE 
"MINUTES" 


ON A LINEAR SCALE 


CALL DRAW(3000- TIME sY¥O9070 pLABEL,TITLD 1030209010909 5e%ei,L AST) 


Y AXIS WILL BE TIME AND LABELLED 


SECOND PEOT 


PLUT 


END OF FIRST PLOT, 


A TIME SERIES REPRESENTATION OF THE ASQB1 


c 


“MAGNITUDE",Y AXIS LABELLED 
CALL DRAW(3 000, TIME sASUs0, O,LABELs TITLB 70309020902 Ue 541, AST) 


EKSUS TIME 


SIGNAL V 
E Y AXIS 


E 


nT: 


LUT OF THE 


BE AP RUe § 
WITH SIGNAL ON THE X AXIS AND TIME ON TH 


THIRD PLOT WILL BE 
& 


OO YQOOOWIUUOU0OY YWUOUOOOOUOUOU YUUOUIO 


CALL DRAW(3000,% TIME sY¥PG,05 OsLABELs TITLE 70402000907 Os dy 4ele LAST) 


fs 


: yap xmas aSiHTSE 
iaiaciisioleeiurassti vat ete ia 


“481 raeaneseye | 


aii 


leieiitalameails 20 ATT o fj ee 


TUSEUG T6LS y2VARA Wi USMITT sia ait UA daudive- 


aia 61 j an i, 4 mn UE ANGT ye QZ d 
as sity Tawe ts uawe lis mar Y Lpeaesaees 
Me JVGTWdan 34 44 Ww a ‘AMS4AAS PINT OD H24. ano Vien 
SIA3E AAA! A. i 3 VIRIAL edad 
S26 J2. G2gN0r1 A A “S21UMi RY O41 2220h2 ba at oo shia <2 bes 3 
t FER ed oO cS Cul 20 «Ged PiteJS2S 4c Oolls OY. an a i eee iS - 
Li a] ¥ Te 4 i ! hig 
Bute 38% Vii ci Aviad ahs ) = ; ® tic { y Pa] =i? 2 rl ‘ 
i> » ao 1 1 
Pat bae 1944 264) - a! em i9He@ ) nA 
WS NA4 ssa AAT ss 4 ie ’ ; 
Bieta yek 
3 * “LA i » Lj iid ( = 
ei wha, = 7 34 
Ot CS’! t | 
i. Cer 
wees € - as ~ s 4 4 


itd 


be." =) 


4¢3¢ J04L ete 


by 


009900300 
OaANOFnor 
aQa0000000 
PEPPPPE 
Qo0000000 
aaaaaaaa 
(aWaWy= Wa eo me A 
2 @ @-o-¢-¢-a~d 


G 
Q 
PLCTTING 


= 
Mud 


FINISHED 


80 


APPENDIX G 
DIGITAL SOFTWARE FOR COMPUTATION OF SYSTEM AMPLITUDE VERSUS FREQUENCY 


BODDDVIDVDSOD A SGACSO OSS ODDDIO VD OOBDDODV DOSOSSCOD00990000009000530 
SD ANP OR DD OAUIFADMDHNOANA Pi ADM DROAN MFA OM DNDN (VFN OF 
UU GN hh ANI NINO AINA FF... 9 OD FY OD SS Se PP PLLA 
OS iS i ES i BS aS I BS IS BS i BS ES a I BS RS i eS SS BS SS OS ES BS BS os Bw ae ES ES RS OS a ES BS eS gO OS 
AWOPVIOOSSSSOSOSSOSOSOSOSSOSNSNSANADDIAOINGSGIVSOSOS9O9NINVOVS9NV9N990095N900909 
See ae a ee a a ee ee ee ee ee eee 


AAAAAAAAAAAA eo SSeS er ee Ee Ee eee Err eer 
PG oe @— oo oO o-oo eo a oe ee oe ee ee oe ng 
(es) 
WwW 
uJ fad 
TZ Uke us 
rALD Bit = 
Ama ~— — 
So) 2 = ~ 
em je Su fo) - 
Za, (Ss) TU = fe) fo) ~ 
—OLiue Oo a if 
Q Dem 9 fos) x, 
An = wd — e Aa 
Wit yore ee > a = 
mnogu act >< oOo wn 
QeNnh a = ~ i cm << 
Za (ea Jee SU) o0U = « o o 
wea DOE eo ow — 6 - 
a— 4 es) oo = i + _ 
ar F~ doe) de) z ga) (Sy ey TE fet ® 
WmemADR NU ~~ oO = .» ~~ oO 
wo 9 => mew Quam WN N fea 
LOQAmbe S| = > -~ wu << = uJ 
WY Us — ms Ge - =~ — 
Q ww m= x eye mm + oN =} 
XY2Z9DFS2Z2 Du “= e& gor +- O C Or ~ 
OmNONH YN - = HQ Qa & aw Ne uw 
— <2 aw Oa NSA a Se) Ww oe 
CAWIWY a5 o270O Qee YW WY = N+ n 
I90=5= 3900 Tea It doe + Ne ” 
>! fer fo af J MND LIQ - = e W — *D — 
Gri a ~w~ww DAY FTN K- - uw +N a. 
RPM Ze a Ker ja) waq-= FOO DV Zz be 
2Z<eaLc ST LUN Sees AO OY i YN ++ =™n 
ay SPS ES Fae am CT4dy = Yy eo = yar 
Fu WwW ne eL< CHM NM YM YY fz Odea 
faa) 2-I2I> z mer “wre oN = = Q Nw mw em TS 
" ft KOOM wu Onrm> Ore See CF le mn 20 0 oO 
7) Zo nN2 3 ooo eset MMN WY Lh rc e+e ONO tue 
7p) meee het ooo eet Foo = w oO ee eNMN UO 
<= ec Oarnw za moo Hee aan = SS — AOMNMDN Bi 
— Ooweye QA NY ~om Lee Coyuy Yn OV xc ~S-NANN Sm 
oO tty Le ty Orww LULL NNW «= NONHN Ro 
- Ae a Be “ NOwe Dee LION QO yO Q -nee% OU 
= trOstOoOunw = DOs ew ey 8 © Ws YZ 1b Nee LS 
220 200 uw aX NU NNN LO <x NeREee YO 
co Ono CORE mare QDee TOOTH oO —~ tere or Ww rh 
- Sworwtwk bt Ome Trea aana ww Wu e— Oe Nw TO 
- ACOKrD COO LwOoOOdtTrTs4o Een O+ OND USO) a} 
ve) On Uy LL. Nn COO FH F®F eLNNNNZ D ew 0 NFEN £ N @D 
o 2) =>- — NOD HO eetetaODnO aw O AORKR+ke O ao wn 
—l =—eoODD ® wWOIM AOU ew ew ON SW LL ~~ +F er UL e ° 
= USER O N= Ow Ae eYeHOnN OO Nr ete - On 
bal Sere w =~ =O OlLaaUx~prpoore YOY YN Fiaewes NW THT 
~- + aiugn2zv2z—ax5 Ud > >< WILL "LAA ” me + aN ws RY NON 
na Cas OS OK OK W<QeenT yj el OW Za Feaen~n oN 2X ONDMD 
POO OWlretewese FAAD eHIINYIMY ZI co ee oC SnD 1 O® Mo Toole Tat] 
=—240 CQO Wa reat Sees A me es es =aza ee ee ae) Q ee 
~~ Cade Wi<t = Sere > Te = seal SS ) riwti~x ON Noow 
K—2 AgqoOnroDArs AANOODODDMDA tt ut as +> PU WW me NP Onp 
Dh EPfOr2DIeY Z2ZSTHw»eseeeieee 48 ZYM N WH ANORN LE oN eww 
OLN Nee wad wWwHIyJJsJIO ew aaaad Le NN ese 
SL BH OROORRRO SBZacsteteaetadee Dae Oral TT te a 
YM Te TMnZewde mw WTO Ww Wo xO bukkake ODP Frnt 
peimes ROOTS INOA AOOCKKKeVeYOO AkK KF O ILEIDMODO OM IOWA 
uw 
Lug >< 0S 
DwOQ 
ae {VE 
nN 


~™N SOU OU OUI UO OVUOOUO YOU OVUUY 


81 


Wocruws 
6 leet 
Cl e094 
pede. oo 
ote & aye 


Oe ewe *} Stim 


spiciacdy aJaeryY QTD ga 


leds hseti, 
eaaremeeetreniihi 
Aete Gedy 8 a 

fb eittipies, 137583 
Te Sicas oe dhe, Brahe 9 hohe WI a Se sg 


Michie GASES » 
ies ate My 


en @-« 


Besse 
eS 


“ 


>= ert 
~ atl de Ht DO 


~- @ bee 


aie She TZ 
7 Plo tH, 


'y 


om 
a 


nA 


+44 yz z 


Fr 
be OD 


LA) Ge Orr e's Sept 


[Sao pee ie @ 
‘ . 


os. - 
é tie 


ary 
Sarvs 


+eSe= 
saa o 
Deo 


os 


! 
1 
E 
j 


Hew Phe 
ie ae 


w 


s LS 
at ot nah 

Sad 

>") 


- 
- 

: > 
ae 
ao be 


$ 
t 
a 
yi 


if 


toes 


99990990990999900000009090900909090990000009009909900 
QDHNOANALPO OM DROANMFO OM DADOANMN FN OEDHDDANMNFPODEDDIDANIY FO 
WVUN0 0000 0000 ORR EERE MHS DD DD DDDDIDDDODHAAAHWARARHODOOSOSD 
PPL PLP PLP PPL FPP LP PT PS PP PT SPE PEA 
9999099990909990099909006009999999999090095699990999909 
AAAAAAAAAAAAAAAARGAAAAAAAAAAAAAAARAAAAAAaAAAaAAaAAAAA 
0200400449 00000004000000000400000402%10040048040000000040, 
qid¢datdqat¢¢@at@ ede dt eqqdt@q@qedt dqeeqdteqqaqdqeqdeteqteqertqetaqi¢ 


e ee ee > 
torr tt + 
SSS Se SN Ok 
= => =e =: 
NNN NN SN 
e tHe te 


DO wh amth te oth th Ott 
OO & FRR RN 
0 ~Newt =m 
He tt HH — Ht He Ht 
OQNDONOWNWL 
OOH + Hie + +t + 
mw es © eof © otf @ 
WON AINE NAN 
~PnwN Rw See 
TON Fe mt He bh be HP 
DN ot Ot +O +L 
Tmet+tO+ Ollie 
aN of + of + ef 
PR en © oN! © ON © 
eee last ain 
Lk oe | | 
NN td mt ed a et 
WL SOO OwiwL oT 


57 


HIon 


FFICIENT FOR THE SELECTABLE 
00% HZ" 


COE 
T O 


1 
MI 


UUOOU 


1¥*G1) 


Q-0——j-2-—]— 
BS i Kine Foe Doe Ke Lonel 
est qGoogow 
= 
ee ed mt ef a 
NOOWLI- 
Oe ed ee ee ee 
ahirniad 
Hun www 
ANH PUVOr 
PPP 
aaaaaada 
ago casas oa 
NNNNMNW 
Po ~ 0 oe Cg 


WITH UPPER FREQ. 


BeeeG hoe FOR SELECTABLE LOW PASS FILTER 


OF O.6 


SSS |S) 


82 


GOL*JJ1) 


*21/4.) 
/ 
/ 


03492 
02779 


e ry a ee ee ee et | 
+ + et 8 OH eM 
SS ~~ NENNE NSM oS 
NI eND eD— = Nw NNN it Ht 6 
DUDMORN N ee Mod oo Te eee Coe Doe Loe 
rs~rrne + +% +O FILHO TS 
OFONNONGH A He DQDOODWLUOOLTHS 
ODDODONDOkKSD KO D+ DW + Uw ww 
Oe eO0 owl wQO + eireoitisytt 


ONODOCOH Hae Hera oO 8 ON oll Ii tl 
SMR OSLO 00-4 1 SI SRO 
eeee eo tN ONO~ww~ww 100 0 
SOQxqeqenFOnW WN aaaa 
HT WE) ME PP sae Fo Kc Coe Lc Lo Lee fe Leen Le Kee oe ee | 
ANMVCHOWU ILOVOWLYO TH AY IUNNNYN 
IOvOWLIOVOWLY De ay ustaad 


eeu sss 


a De 
We “O44R 


ace Ay 


rertin Sects 38 Bit ABA one A iT 


oVGRS B99) WITS BHL44 COSR a) Piwata= 


Euan 


74 eris}) a4 73 


ool EE co 
Bee de ee 


ie | 2U = roe wile =) 


A & at Los cay 
e® Lo4 iat reise 


: 


=) 


ae) as ae 
4 


SDd0OCoS DON VOSS NOVO SOON SDD VV IA VN NVONCSOSINIVI OVO ODOC ON COO 30090900 
DM DDO SQN FPN OM DO DOSNUMFD OM DADANMFHOMDDOAN MOMSEN OM DRHOANM 
DOD eS hh ANI NO IN ASO IN 9 FF COVED EOVETY OO CH OF) SPS SS SPLAT 
a ral alta Via lla lla lla lialialvalvalvallallalalralalalrabalralial alralalialalalalialalalalalial al alepalaliallalelalatalay 
oso og BDO900C COO 9N 9000 9000 9089N90000 0000900090090 0000909090909 0 
AWaAGAaAAAARAAAAAAARAAAAAMaAARAAAAAAAAAAAAAAAAAaAaAAaAAAaAA 
AQAAAAAALAAAAAAAAAAAAAAAAAAAAARAAAAAAAAAAAAAAAAAaaAA 
Tot it oe I gf I fe SI ef es I I I et It td Qe ee et e@ taster 


GG1l*JJl) 


<--> 
=< 4 
oe Ler) 
234 
+#O%* 
a0 
OO 
On 
—~—w\ w= 
—~ON-4 
q+ 5 
> To soe Lae BP) 
% QU 
“AOWWR 
See eo 
Ons 
—ndoUO ~~ 
S30O™N 
_—Fe +e em 
em Lente Lae ame 
2) te UL Le 
A uwWiLL UW LL 
te st it ot 
4 rd 
i 2 © doale) 
o> Go oe daale) 
~we et gg ee ee 
wut wy 
SANS 
0 0.0.0.0 
aAaAaaA oO 
Jaye 
NMNWY) 
anonnao 


STORAGE REGISTERS 


INITIALIZE 


FINISHED COMPUTING CGEFFIECIENTS, 


EQUATIONS FOR SIMULATED SIGNAL 


>) 
jen) 
j— 
— 
ud 
“ 
Pe) 
aw 
(j@) uw 
>) = 
ie) Ta) 
My — 
» oO 
= uw 
1 id ees (a 
N39 2eOOOwW 
wonnhw>D ww 
OmiwnremZz 0 
OND=95G7— << 
Qa ~~ mespes 
2~O0O9G0O2Z O 
OxeOAMs=O — 
QR S$ eNO YN 
fo) 
(>) 
NI 
OUUY 


—— 

ak 

+ % 

ee 

NN 

+ + 

ej 

QQ 

Ly LL 

WLU fo) 

NY oo 
2QO Oooo0o00000 NNO 
NO wee AQUOS VO POD DRI TYP UNO OOO ath 
MMNNDOOCDOCO00CSO £8 2e & @ &@ & 2 & © ot & 
Se ee ee ee lelel=lelelelelerlelin | 
ANDOOOOOCOCOCOOH HH UH UW nb epee 

DD WW EE TD) ee mm ee ce oe a oe ee a mm rt pg 


DORN NDO HO Dee St Rt St FNIOORNIM PUN 
CN me me rr re re re et ee ee er ee ee OO} fl} ltl lt 
WeANIWIIVGIISIBSIBIBUIBBIBISBGIIIN AS Ht 
Cpe me LLP Ly Tn) i | 
ALT LLCLCLCSC LLL CA EK LYK CLOORNM SN 
fam Ya Wo We ye SUE ee ee ee a anyon) a] 


83 


POODEL OOS Ow1o 


hate wae e ; oy 


“YT pose ae rvs 


Ne ene et) ere 
¢ a BO de deb fri 
SPret srt seg 


: 


{} 


ton 


gc 
oe 


o°  EOeh Oe? ber THMVEED 2peeT 
HEC COLL 1? \CORELIECIEWL?? JuELtwrixe 


oe 
£8 
» @ 
a | = 
> B&B 
~<- 
j tf 
ae ? 
ae <— 
‘ <; 
_— —— (=~) e— 5a A “Mure ot, o 
ee > ee iene an. ® 
- on 24 & ©» ace 6 Laine { 
a, 1 St Cease, S 
> ee ee i a re —- 4 


lm ie oe — o% 


ny 
bi & 
re. ies Fe, A ‘ee roc Carin 
: (Matyas ea rages i 


‘ ae Cae ~~ eek a 


a plas 


eager 
Sa7 te 


oe 


BSDDOOCOP OSGOOD OO DO VOSS SSO VO Sg 9oV a gg 0eq0 00 90090 909CNOVOSO9NNO 
FIDO DDO SQUID Oh DRO SAIM FINDER DAO ANOLON OHS DHOANUOASNOMODROA 
LAV LIVLAVLAA LILA SO 0 WO 00 00.90 ORE Re Re Reh mM RM MH DDO DO DHDDDDAAVAADADIARDADAOO 
WoValtalralralral alalralalraltalalralalal al alrallial allel alralalialialalalialalalalialialalial al alialallalal alaliabele) 
Qcoo0o0qo0 0090000 009NDIGDOV OV OV CO NVN0N9NSCON SBS SEN S98909gN9N909N990 
AaAGWaAAAAAAAAAAAAAAAAAaAAARAAAgRAgRaAaAaAAAAaAAaAaAAaAaAnaAAaAaAAGa 


AAAAAAAAAAAAAGBAAAAAAAAAAAAAAAAAAAAAAAAAAAAAARAAAe 
{<4 tet ft ef OI SI ef ef I tt I I< 
Zz uw 
— Qa 
Ze = = 
ud a ¢ ~ a = 
Aime - — a 
230 —_ = = 
qn" = ie |] = 
Oo ~e 4 n () 
wmnpnuI w em a 
ome Ky Wa ~ + be) 
ago = So + ae 
Revue » o — 
WMNULoet—— + a 
ORYMImie st =) Le -_ 
a= Io a ” o>) Ww 
QwZOornmM FT a e 
Wb AM Wh + “n fo) 
= Deiiet = ee) = 
ZyL2Z=a + o— = uw 
WwW Le = =~ =— mn if 
SanzoS N a N ~—— <4 = 
=arF=zOoOoO ~ ~~ o~ > = 
Oo <0O>-—-~ ~~ AQ HT} 
SAS) 25 — ta uy 
— < — + > CD Fe) 
New W + — = M+# <= ”n 
—Ie=r uw + mo -«S NN *ON [a of es) 
SRS ON rite =~ avo Ww =) 
=S=Oz ua A aA nw ew IC > 
= tw OFTu £ zt > YS =z a 
- “LxZzZ2ze~ uL nN awx< <n a 
— S4Om i a rt CN+ +0 iu 
OLN>oOr + + NRE =+ = 
Led ec sOZKe = -] dep N= ad = 
=z QwtSe ol N +3Q 3 = 
a | Tawseaw = — —s—Tr He YNZ N 
4 CReaAwvwooao~ ~~) - Own wo a 
Uy MiSs >) 3 eK THN Or ae tay WO 
5 rPE=OKyYII Sawer Se Aa 
= 4IaZUDNOY wt mews Ht OO af Ow 
4 — oe Lee Ot at mt ad > OS R=) nS 
=a LL. Ou ~ >< Stet 0 4 Sz -- —DOe 
~ Ow eA a Ot <0 OF Qu ANN fom | an) 
ie | weoowWd = XT ew > ear — te — oi] 
Ls e CTlaadt2Z2ere w AN See< awn so Led +e +f = 
ecoET 2M faa) tite 1 tHa y -— ~~ — 
uw Oo Ug Ly uy + +a <0" +0 — Ome fe Ae) 
+ Vole Ueae u0 = mn Tt FOIm + =O (o> he bee |} ~~ > 
x ~ OFQOuNISt =) *O sr@= Oa Oww FKOOW 
(o Reet fw wwe IOI wer et + uJ MONOID ABawDdD 
Satter UU >) Oom-+ 1n>—~— NO Or NODA 
WPuUUeeN ~~ GOMwWReRw «= ~~ em IK IO NN ~ Ome TOWW 
ANAPING ome Lf LEDOQ N= % pee eH o Ow~w NI 
ies Oat fea OO Ne Hewes Ons ww aun Ore NORL 
me OOOO Ott mt mt tes te 3 OO + he NOK GHtO= 
Seat teat ZS QL weer we et te HOS x Ja 4 uND N20 @ 
eeoeroeseommvbl FEY NZO LOOT >€4HS>>00 eo Ow OOTD=Y TOrx< 


RRR REeROe tt Neos LL Sort I MIKO HO 
Saag 4dndyS = Oe ES OHH MH I I 


j= ee QO ZWD eM 
+IZz O00 MONnDMNIDO~ 


Oe ee 0 8 emmy SOQ 2S NH am 1 -NNO= 00 When WtZNnwe 
LANOSINO—= Res LL DUWWYE Oe eNweree I oOaeH wb BCROCr-OMw<«t 
SF ram rad Dt en Pe mm se? es LJ NILEO SLA ewrwr~ i it im NANNGRH-OK-= 


SOLLLL LL LLILiLOfes CPSU Of ORR Ree>Saaaow 
tt et eget ee I ke kR— ete, | +} > > OO KY YH 
Sed * 


ZW TPNZRIHe 
r2z DzODTFOstAOLO 
Het NNONNOCOSYL 


—} 
oO 
mo 


100 CONTINUE 
E 
D 


4001 


QUOOUOUOUUYU OO YUOUO 


84 


: 
a 


— i 


ae 
ae git, 

Celene. 
er ee) eT, 

ty ¢ ¢ 

eC oR 
hig, 
BSCS 


are 


ever? 


STIS bOVEUETeEOITT 


hem i 
= od 
‘ «a 
—_ 
= = ow ld. 
x —— i 
owe %, ~- 4 
i, ~ ~~ 
é i tee 
4 xt — 
a — - - 
- = ww 
e 


a - ff is - =) 
—_— = som > 
\ fi ° . i Ce j - 
—£ a SS - 
‘ i=T > bd i oe Sel = 
“ or, + * i; Coaety -— 
a a - ~—Slee Sp -~ 
. “ ry on SD » 
~ - ~ it a ite TS = 
ae en we 23x Pee t 
~ > © ose a i © > ( P< 
a i Oe - = New ¥ 
we or | Cite... = 
“i nll ey oh! be = 
A, Maciel it oy : i = 
Li aid 
ee Ae i é dav | - 
ae <i me 2 qLit o Corea: <c 
tb et a ar re =) sa ‘© tai 
a fi fa 4 « = — im 
“4 ioe — [a ~ 
i * ath SS 
Ai, , f he a a «#- 
ey i P i TS eo 
= ; aa — _ : ae im . 
: * ’ £ 
tity me g 5S ww 
des, gem + 4 ~~ 


« - ‘ =< } os ae a 
om = = > pet ea 
Py * 1 po and fo. Ure Pio 
be ith git “a. . é Sc Ow f mete og 


: re ce eee hee 
a. th om -_ ‘ wrk ‘ i ee ee 


2 -— = te = 
—e 
Ie ars . 
i 1 ten P * te t ~ Iw 
" 
iw P. ‘i 7 - 
i pe > € joe 
-_ 
Ps ay 4% 
4 


e0ooosos 
Narn oreo 
as0900S530 
0 0.0.0.0 .0 0 
(ole le lolale lo) 
aaaaaaa 
aA &®A0,00, 0, 
taqaatad 


FINISHED PLCTTING 


85 
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